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Nonequilibrium molecular dynamics simulations often use 
mechanisms called thermostats to regulate the temperature. 
A Hamiltonian is presented for the case of the isoenergetic 
(constant internal energy) thermostat corresponding to a tun- 
able isokinetic (constant kinetic energy) thermostat, for which 
a Hamiltonian has recently been given. 

Thermostats are modifications to the equations of mo- 
tion of a classical system to simulate thermal interac- 
tion of a system with the environment. The Nose-Hoover 
thermostat is used to simulate fluctuations in energy of 
an equilibrium system corresponding to the canonical 
ensemble of statistical mechanics, and the Nose-Hoover 
and Gaussian thermostats, among others, are used to re- 
move heat from a system driven by external forces into 
a nonequilibrium stationary state ||l|. There has been 
a recent interest in thermostatted equations of motion, 
focussed on the symplectic structure of the equations of 
motion, and the related pairing of the Lyapunov expo- 
nents. Both a Hamiltonian and pairing of Lyapunov ex- 
ponents are known for Nose-Hoover and Gaussian isoki- 
netic (GIK: constant kinetic energy) thermostats [0-0. 
Numerical evidence against pairing (and hence the ex- 
istence of a Hamiltonian) are discussed in M for the 
GIK thermostat applied to shearing systems and in g 
for the Gaussian isoenergetic (GIE: constant internal en- 
ergy) thermostat. The latter paper does, however show 
that in a special case of the GIE thermostat, involving 
one rather than two arbitrary potentials, the Lyapunov 
exponents are paired. The purpose of this short note is 
present a Hamiltonian for this case. 

The GIE thermostat has equations of motion of the 
form 
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where $('^^*' is the external driving potential, <I>(™*-' the 
interparticle potentials, and a is the thermostat term 
which ensures that the equations conserve internal en- 
ergy E = X]iPf/(2TOi) -I- <I)(™*). The equations reduce 
to no thermostat when ^t'^^*) = and to GIK when 
(j)(mt) _ Q^ ^ more general example of a limit involv- 
ing only one arbitrary potential is the case $('^^*) = 7<I>, 
<j)(mt) _ j^j^ _ ^)(j)^ leading to the equations 
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which conserve energy E = X]iPf/(2?^i) + (1 ^ l)^- 
Here, 7 effectively controls the strength of the thermostat 
from no thermostat (7 = 0), to the GIK thermostat (7 — 
1). For any 7 the Lyapunov exponents are paired [^, 
suggesting the existence of a Hamiltonian. 

Following the GIK case , the conservation law is en- 
forced by setting the numerical value of the Hamiltonian 
equal to the conserved energy, assigned the value zero by 
a shift in the potential energy. This allows the kinetic 
energy term in the denominator of (II]) to be replaced by 
minus the potential energy (note <I> < 0), 
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Another aspect of a Hamiltonian description of ther- 
mostatted systems is that in the physical variables (x, p) 
there is a phase space volume contraction rate propor- 
tional to a, while in the canonical variables (x, tt) phase 
space volume is conserved. This means that tt must 
be greater than p by a factor equal to exp (J adt^ = 
$|7/(2(i-7))^ Multiplying the zero energy by an arbi- 
trary power of |<I>| we have 
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which, combined with the constraint Hp = and 
the identifications dt = |$|-t/(2(i-7))+/3(;;\ and p, = 
|$|-7/(2(i-7))7r, leads to the equations of motion, (|). 
Interesting cases are /3 = 7/(2(1 — 7)) for which there 
is no time scaling, /? = has a certain simplicity, 
/3 = —7/(1 — 7) yields the familiar form of kinetic plus 
potential energy, and /? = — 1 for which the Hamiltonian 
is that of a geodesic in a conformally flat space, see [g|. 
The author is grateful for discussions with W. G. 
Hoover. 
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